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Abstract 

Discrete flavour groups have been studied in connection with special patterns of 
neutrino mixing suggested by the data, such as Tri-Bimaximal mixing (groups A4, 
S4. . .) or Bi-Maximal mixing (group S4. . .) etc. We review the predictions for sin0i3 
in a number of these models and confront them with the experimental measurements. 
We compare the performances of the different classes of models in this respect. We 
then consider, in a supersymmetric framework, the important implications of these 
flavour symmetries on lepton flavour violating processes, like fi ^ and similar 
processes. We discuss how the existing limits constrain these models, once their pa- 
rameters are adjusted so as to optimize the agreement with the measured values of 
the mixing angles. In the simplified CMSSM context, adopted here just for indica- 
tive purposes, the small tan/3 range and heavy SUSY mass scales are favoured by 
lepton flavour violating processes, which makes it even more difficult to reproduce 
the reported muon g — 2 discrepancy. 
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1 Introduction 

Neutrino mixing [1-6] is important because it could in principle provide new clues for 
the understanding of the flavour problem. Even more so since the neutrino mixing angles 
show a pattern that is completely different than that of quark mixing. The bulk of the 
data on neutrino oscillations are well described in terms of three active neutrinos. By now 
all three mixing angles have been measured, although with diflFerent levels of accuracy (see 
Table 1 [7]; see also Ref. [8]). In particular, we now have flrm experimental evidence for 
a non-vanishing value of the smallest angle ^13 and a rather precise determination of its 
range (see Table 2 [9-12]). 

Models of neutrino mixing based on discrete flavour groups have received a lot of at- 
tention in recent years [14-19]. There are a number of special mixing patterns that have 
been studied in that context. The corresponding mixing matrices all have sin^^23 = 1/2^ 
sin^ ^13 = 0, values that are good approximations to the data, and diflFer by the value of 
the solar angle sin^0i2. The observed sin^^i2, the best measured mixing angle, is very 
close, from below, to the so called Tri-Bimaximal (TB) value [20-24] of sin^ ^12 = 1/3 
(see Fig. 1). Alternatively, it is also very close, from above, to the Golden Ratio (GR) 
value [25-28] ^m^ 9^2 = ^ = ^ ^ 0.276, where = (1 + y5)/2 is the GR (for a 
different connection to the GR, see Refs. [29,30]). On a diflFerent perspective, one has also 
considered models with Bi-Maximal (BM) mixing, where sin^6^i2 = 1/2, i.e. also maximal, 
as the neutrino mixing matrix before diagonalization of charged leptons. This is in line 
with the well-known empirical observation that 612 -\- 6c ^ 7r/4, where 6c is the Cabibbo 
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Table 1: Fits to neutrino oscillation data from Ref. [7]. The results for both the normal 
and the inverse (in the brackets) hierarchies are shown. 



Quantity 


sin^ 26'i3 


sin^ 613 


T2K [9] 
MINOS [10] 
DC [11] 
DYB [12] 
RENO [13] 


0.111°- (0.141°:-) 

0.0411°:°^! (0.0791°:°^^) 

0.086 ± 0.041 ± 0.030 
0.092 ± 0.016 ± 0.005 
0.113 ±0.013 ±0.019 


omst^oll (0.0361°:°^^) 
omot'Zl (0.0201°:°}^) 

0.0221°:°}^ 
0.024 ± 0.005 
0.029 ± 0.006 



Table 2: The reactor angle measurements from the recent experiments T2K [9], MINOS 
[10], DOUBLE CHOOZ [11], Daya Bay [12] and RENO [13], for the normal (inverse) 
hierarchy. 



angle, a relation known as quark-lepton complementarity [31-51]. Probably the exact com- 
plementarity relation becomes more plausible if replaced by 612 + 0{9c) ^ 7r/4 (which we 
call "weak" complementarity). One can think of models where a suitable symmetry en- 
forces BM mixing in the neutrino sector at leading order (LO) and the necessary, rather 
large, corrective terms to 9^ arise from the diagonalization of the charged lepton mass 
matrices [31-57]. Thus, if one or the other of these coincidences is taken seriously, models 
where TB or GR or BM mixing is naturally predicted are a good first approximation. 

Exp 

I 



GR TB BM 

Figure 1: The values of sin^ Ou for TB or GR or BM mixing are compared with the data 
at la. 

In the following we will mainly refer to TB or BM mixing which are the most studied first 
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approximations to the data. The simplest symmetry that, at LO, leads to TB is A4 while 
BM can be obtained from 5^4. A4 models have been studied widely (for a review and a list 
of references, see Ref. [14]). At LO the typical A4 model leads to exact TB mixing. The LO 
approximation is then corrected by non-leading effects. Given the set of flavour symmetries 
and having specifled the fleld content, the non-leading corrections to TB mixing, arising 
from higher dimensional effective operators, can be evaluated in a well-deflned expansion. 
In the absence of speciflc dynamical tricks, in a generic model all three mixing angles receive 
corrections of the same order of magnitude. Since the experimentally allowed departures of 
9i2 from the TB value, sin^0i2 = 1/3, are small, numerically not larger than 0{Xq) where 
Xc = sin^c, it follows that both 9is and the deviation of 623 from the maximal value are 
also expected to be typically of the same general size. The same qualitative conclusion also 
applies to models with GR mixing. This generic prediction of a small ^13, numerically 
of (9(A^), is now confronted with the most recent data. The central value sin^ia ~ 0.15, 
from Table 1, is between (9(A^) - (9(0.05) and 0{Xc) ^ 0(0.23). Since Ac is not that 
small, this gap is not too large and one can argue that models based on TB (or GR) mixing 
are still viable with preference for the lower side of the experimental range. 

Of course, one can introduce additional theoretical input to improve the value of 9is (for 
an updated list of recent models of this kind see Ref. [58] and references therein; for recent 
models in the speciflc context of holographic composite Higgs see Refs. [59,60]). In the 
case of ^44, one particularly interesting example is provided by the Lin model [61] (see also 
Ref. [62]), formulated before the T2K, MINOS, DOUBLE CHOOZ, Daya Bay and RENO 
results. In the Lin model the A4 symmetry breaking is arranged, by suitable additional 
Zn parities, in a way that the corrections to the charged lepton and the neutrino sectors 
are kept separated not only at LO but also at next-to-leading order (NLO). This way, the 
contribution to neutrino mixing from the diagonalization of the charged leptons can be of 
(9(A^), while those in the neutrino sector of 0{Xc). In addition, in the Lin model these 
large corrections do not affect 612 and satisfy the relation sin^ 623 = l/2 + l/v/2cos5 sin^ia, 
with S being the CKM-like CP violating phase of the lepton sector. Thus, in the Lin model 
the NLO corrections to the solar angle 6^12 and to the reactor angle 613 can be naturally of 
different orders. 

Alternatively, one can think of models where, because of a suitable symmetry, BM 
mixing holds in the neutrino sector at LO and the corrective terms for ^12, which in 
this case are necessarily rather large, arise from the diagonalization of charged lepton 
masses [31-57]. These terms from the charged lepton sector, numerically of order (9(Ac), 
would then generically also affect 613. The resulting value could well be compatible with 
the present experimental values of 613. An explicit model of this type based on the group 
S4 has been developed in Ref. [52] (see also Refs. [53-55]). An important feature of this 
model is that only 612 and 613 are corrected by terms of 0{Xc) while 623 is unchanged at 
this order. This model is compatible with present data and clearly prefers the upper range 
of the present experimental result for 613. 

In this work we discuss three possible classes of models: 1) typical A4 models where 9 13 
is generically expected to be small, of the order of the observed departures of 9^ from the 
TB value, and thus with preference for the lower end of the allowed experimental range. 
2) special A4 models, like the Lin model, where 9i3 is disentangled from the deviation 
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of sin^0i2 from the TB value 1/3 and can be as large as the upper end of the allowed 
experimental range. We give a general characterization of these special A4 models where 
the dominant corrections to TB mixing do not arise from the charged lepton sector but 
from the neutrino sector. 3) Models where BM mixing holds in the neutrino sector and 
large corrections to 612 and ^13 arise from the diagonalization of charged leptons. The 
value of 9is could naturally be close to the present experimental range. In each of these 
possible models the dominant corrections to the LO mixing pattern involve a number of 
parameters of the same order of magnitude, ^. We discuss the success rate corresponding to 
the optimal value of ^ for each model, obtained by scanning the parameter space according 
to a similar procedure for all three cases. We argue that, while the absolute values of the 
success rates depend on the scanning assumptions, their relative values in the three classes 
of models, provide a reliable criterium for comparison. We find that, for reproducing the 
mixing angles, the Lin type models have the best performance, as expected, followed by 
the typical A4 models while the BM mixing models lead to an inferior score, as they most 
often fail to reproduce 612. 

We then discuss the implications for lepton fiavour-violating (LFV) processes of the 
above three classes of possibilities, assuming a supersymmetric context, with or without 
See-Saw. The present bounds pose severe constraints on the parameter space of the models 
(for a recent general analysis on model-independent fiavour violating eflFects in the context 
of fiavour models, see Ref. [63]). In particular, we refer to the recent improved MEG 
result [64] on the /x ^ 67 branching ratio, 5r(/x 67) < 2.4 x 10"-^^ at 95% C.L. and 
to other similar processes like r ^ (e or /a)^. One expects that lepton fiavour-violating 
processes may also have a large discriminating power in assessing the relative merits of the 
different models. We have studied this by adopting the simple CMSSM framework. While 
this overconstrained version of supersymmetry is rather marginal after the results of the 
LHC searches, more so if the Higgs mass really is around uih = 125 GeV, we still believe 
it can be used here for our indicative purposes. We find that the most constrained versions 
are the models with BM mixing at LO where relatively large corrections directly appear 
in the ofi[-diagonal terms of the charged lepton mass matrix. The A4 models turn out to 
be the best suited to satisfy the experimental bounds, as the non-diagonal charged lepton 
matrix elements needed to reproduce the mixing angles are quite smaller. An intermediate 
score is achieved by the models of the Lin type, where the main corrections to the mixing 
angles arise from the neutrino sector. Overall the A4 models emerge well from our analysis 
and in particular those of the Lin type perhaps appear as the most realistic approach to 
the data among the discrete fiavour group models that we have studied. As for the regions 
of the CMSSM parameter space that are indicated by our analysis, the preference is for 
small tan /3 and large SUSY masses (at least one out of mo and mi/2 must be above 1 TeV). 
As a consequence it appears impossible, at least within the CMSSM model, to satisfy the 
MEG bound and simultaneously reproduce the muon g — 2 discrepancy. 

The paper is organized as follows. In Sect. 2 we discuss the models and their predictions 
for 9is to then compare them with the data. In Sect. 3 we list the efiFective operators that 
induce lepton fiavour violation and derive their contributions to the measured quantities. 
We then apply the general formalism to the specific models and observables. In Sect. 4 we 
derive our conclusions. 
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2 Models 



We consider models invariant under a flavour symmetry group Gf. At the LO the 
lepton mixing arises from a mismatch between the residual symmetries Ge and G^, of 
charged lepton and neutrino sectors, respectively. In this approximation charged leptons 
and neutrinos acquire mass from two independent sets of flavons, $e and whose VEVs 
preserve two Abelian groups: Ge = {n > 3) and Gjy = Z2 x The groups Ge and Gj^ 
can be subgroups of Gf or, as a result of the speciflc fleld content of the model and of the 
LO approximation, they can contain some accidental symmetry and generate a group G 
different from (and possibly larger than) Gf} Lepton mass matrices can be expanded in 
inverse powers of the cut-off scale A 

rrie mf^ + 5m^^^ + . . . {2 1) 

rrijy m^^^ + 5m)}^ + . . . 

where 5m^J^l/m^^l = (9(($e,z/)/A) and dots stand for higher-order terms. The LO contribu- 
tions TTZe^^and m[/°^ depend only on ($e) and respectively. They are invariant under 
Ge and G^: 

Here gei and gj^i are the elements of Ge and Gjy and p denotes an irreducible three- 
dimensional unitary representation of the group G generated by Ge and G^.^ Both Ge 
and Gjy are Abelian and the matrices p{gei) and p{giyi) can be diagonalized by two inde- 
pendent unitary transformations Qe and 

P{9m)diag = ^1 p{9ui) , p{gei)diag = p{9ei) , (2.3) 

and the mixing matrix Upmns reflects the misalignment between the two bases: 

UPMNS = ^l^u . (2.4) 

It is immediate that the mixing matrix is independent of the choice of the basis. We 
choose to work in the basis where p{gei) and me^^^me^^ are diagonal: = L Beyond 
LO, the prediction for the mixing in eq. (2.4) is modifled. In general the VEVs of $e and 
are corrected by terms of relative order {^e,h)/^ and do not preserve Ge and Gi, any 
more. Moreover higher order operators contribute to lepton masses without respecting 
the LO residual symmetries. The NLO corrections are suppressed with respect to the LO 
contributions by the ratio between the flavon VEVs ((^e), {^y)) and A. Depending on the 
agreement of the LO approximation to the data, ($e)/A and ($^)/A will typically range 
between and A^. 



^ In the present paper we concentrate only on the TB, GR or BM patterns, originated by the mismatch 
between Ge and Gy. However, interesting deformations of these patterns, all predicting already at the LO a 
non- vanishing ^13, arise considering the finite modular groups Fat, > 1: in Refs. [65,66], a comprehensive 
analysis is presented for an arbitrary Gg and G^^ = Z2 x Z2. In Refs. [67-69], a more general study has been 
presented, where the residual symmetry in the neutrino sector is Gy = Z2, while the other Z2 component 
arises accidentally. 

^Later on we will include in m^^ also higher-order contributions satisfying the property (2.2). 
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2.1 A4 Models 

We refer to SUSY models based on the flavour symmetry Gf = A4 x Gaux^ where 
the Gaux factor depends on the speciflc reahzation [61,70-72].^ The group A4 can be 
generated by two elements S and T satisfying 

= {STf = = 1 . (2.5) 

The irreducible representations of A4 are a triplet 3 and three inequivalent singlets 1, V 
and V\ In the triplet representation S and T can be chosen as: 



T= , 5=- 2 -1 2 , (2.6) 

\00cj/ ^\22-l/ 

where uj = e^^^/^. Under A4 the electroweak SU(2) lepton doublets / transform as a triplet, 
while the electroweak singlets e^, fi^ and and the electroweak Higgs doublets as 
singlets. In the flavon sector both $e and always include a triplet, but they can also 
include additional singlets. At the LO and in the exact SUSY limit the VEVs of $e and 
are determined by two separate sets of equations and satisfy at LO: 

r'($e) = ($e) , S{^,) = . (2.7) 

The transformation can coincide with the T generator of A4 [70,71], or can represent 
an accidental symmetry of the charged lepton Lagrangian still satisfying T'^ = 1, as in the 
models of Refs. [61, 72], where = ujT. The charged lepton mass matrix rUe is given by 

rUe = m(°) + Sm^^^ + . . . (2.8) 

where 

mf=t-J \v- (2.9) 




y- 



T 



Here va is the VEV of i/^, yf (/ = e^fi^r) are dimensionless quantities and is a small 
parameter that breaks the flavour symmetry A4. At the LO the charged lepton mass 
matrix m^e \ only depending on (^e), is diagonal and invariant under the action of the 
transformation T': 

T'^mf t^(o)T' = t^(o) . (2.10) 

We have = Z3, generated by T^ The hierarchical pattern ?/e <C <C i/t can be 
reproduced by requiring that operators of increasing dimension contribute to and ^g. 



^The present analysis applies also to models based on flavour symmetries that contain the group A4, 
such as 5*4: few examples can be found in Refs. [73-76], where the TB pattern is predicted at the LO, while 
at the NLO the mixing is corrected in a similar way as we are going to discuss in this section. Moreover, 
our analysis applies also to the model based on the flavour group T^, described in Ref. [77], even if does 
not contain A4 as a subgroup. Although we could generalize our analysis to a brother class of symmetries, 
we focus on the flavour group A4 that represents the minimal choice in terms of dimensions of a group. 
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respectively ^. This can be achieved either by means of a Froggatt-Nielsen U{1) symmetry 
contained in Gaux [70,71] or through some discrete components of Gaux as in Refs. [61,72]. 

This class of models can be realized both with and without a See-Saw mechanism. In 
the first case there are three right-handed neutrinos transforming as a triplet of A^^ while 
in the second case the source of neutrino masses is a set of higher dimensional operators 
violating the total lepton number. In either case the light neutrino mass matrix is given 
by: 

= m^^^ + Sm';}'^ + . . . (2.11) 

where 

I X y y \ 

m[°) = y x + z y- z . (2.12) 
\y y- z x + z J 

The parameters x, y and z are quadratic in the VEV of H^, and inversely proportional 
to the scale associated with the violation of the total lepton number. The LO term mi°^ is 
invariant under the Z2 x Z2 symmetry generated by S, eq. (2.7), and by 

/I 0\ 

A23 = 1 . (2.13) 
\0 1 oj 

Indeed, 

S^m<;^^S = mf and Alm<^^^ A23 = m^^K (2.14) 

The Z2 symmetry generated by the matrix ^423 is an accidental symmetry. The matrix 
m^/^^ of eq. (2.12) is the most general one invariant under both S and A23. In the minimal 
formulation of Refs. [61,70,71] the parameters x, y and z are not independent. If neutrino 
masses are generated via the See-Saw mechanism [61,71] they are related by 

z=^-^^. (2.15) 
Ay - X 

If neutrino masses are parametrized directly through a higher dimensional operator [70,71] 
we have 

z = -{x + 2y). (2.16) 
At the LO me^^ is diagonal while m[.°^ is diagonalized by 



V 3 V 3 ^ 



TB 



1 /i _ /I 
6 V 3 V 2 



V -a/- a/- a/- / 
\ Veya V2/ 



(2.17) 



This holds for any value of the parameters x, y and z since Utb simultaneously diagonalizes 
both S and 7423- 



^Here, we include in mi^^ all contributions arising from the LO (^e), independently of the dimensionality 
of the corresponding operator. 
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The departure from the LO approximation depends on the subleading contributions 
Sm^e"^ and 5ml^\ which can vary for different models. In aU models considered here [61,70- 
72] the NLO correction to the charged lepton mass matrix is not invariant under and is 
of the following type: 

/ 0{ye) 0{ye) Oive) \ 
5m^^^=vA 0{y,) 0{y,) 0{y,) U^, (2.18) 
V 0{yr) 0{yr) 0{yr) ) 

where is a small adimensional parameter given by the ratio between a flavon VEVs and 
A. The transformation to diagonalize mg is V^rrieUe = mf^^ with 

Ue= -cT2i 1 cisH (2.19) 

-cTsC 1 J 

where c^2 5 ^13 C23 are complex parameters of order one in absolute value. We discuss 
the NLO contribution to by distinguishing two cases. 

2.1.1 Typical A4 Models 

In "typical" A4 models [71,72], the NLO contribution Sm^J"^ is a generic symmetric 
matrix with entries suppressed, compared to the corresponding entries in m^u \ by a relative 
factor (^', of the order of the ratio between a flavon VEV and A. This occurs both with 
and without the See-Saw mechanism. The generic transformation that diagonalize rriiy is 
Utb ' where 

U,= -dili' 1 c^3^' , (2.20) 
K-clli' -dili' 1 I 
where 0^2? ^13 C23 are complex parameters of order one in absolute value. Barring a 
flne tuning of the Lagrangian parameters, in these models the suppression factors ^ and 
^' are expected to be of the same order of magnitude. For example, beyond the LO the 
equations satisfled by ($e) and ($^^) are no longer decoupled and the corrections to the 
LO flavon VEVs are of the same size, for both $e and All the elements of the mixing 
matrix get corrections of the same size C ^ C'- We expect^: 

sin^ ^23 = ^ + 7^e(c^3) i + {Ueic'i^ - ^2 7^e(c^3)) ^ 

sin^ e,, = \- \ne{cl, + cl,) i + ^ ne{di,) ^ (2.22) 



sin 9 



13 



3^ (c^2 - c^3) + 2^3 (V^di^ + c^g) 



^Eq. (2.22) is a particular case of the general parametrization presented Ref. [78]: 

sinl923 = -^(1 + a), sinl9i2 = -^(l + s), sinl9i3 = , (2.21) 

with a, s and r real numbers. The expressions in Eq. (2.22) show explicitly the dependence of the NLO 
mixing angles on the corrections from both the neutrino and the charged lepton sectors. 
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We see that to reach the central value for the reactor angle in agreement with the value 
reported in Table 1, the parameter ^ should be of (9(0.1). A precise value can be found 
studying the success rate to reproduce all the three mixing angles inside the corresponding 
3a ranges, depending on the value of ^. As shown in Fig. 2, the value of ^ that maximizes 
the success rate is 0.076(0.077) for NH (IH). The corresponding value is ^ 8.5%, which is 
not large but not hopelessly small either. 



Figure 2: Typical A4 Models. Success Rates as a function of the parameter ^. The qj^ 
parameters that multiply ^ are treated as random complex numbers with absolute values 
following a Gaussian distribution around 1 with variance 0.5. In Cyan the NH and in 
Orange the IH. The value of ^ that maximizes the success rate is 0.076(0.077) for NH (IH). 

We analyze quantitatively the expressions in eq. (2.22) and their correlations in Fig. 3: 
in the plots on the left (right), we show the correlation between sin^ ^13 and sin^ ^12 
(sin^^23)- The parameter ^ is taken equal to 0.076. The qj^ parameters that multiply 
^ are treated as random complex numbers with absolute values following a Gaussian dis- 
tribution around 1 with variance 0.5. In the plots we show only the NH case. The IH case 
is similar. 

As we can see, the plots are representing the general behaviour of this class of models: 
sin^^i3 increases with but correspondingly also the deviation of sin^ ^12 from 1/3 does. 
As a result, even for the value of ^ that maximizes the success rate, the requirement for 
having a reactor angle inside its 3a error range corresponds to a prediction for the solar 
angle that spans all the 3a experimental error bar and is often not even in agreement with 
the data. 

2.1.2 Special A4 Models 

In these models the accidental symmetry A23 of the neutrino sector is broken by a 
relatively large amount so that, in first approximation, the residual symmetries of the 
charged lepton and neutrino sectors are those generated by and respectively. At the 
LO and in the chosen basis m\me is diagonal while is invariant under S\ 




OJM 0.10 0.S5 020 0.25 



S'^rrij.S = m- 



(2.23) 
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(a) Correlation between sin^ ^12 and sin^6>i3. (b) Correlation between sin^ ^23 and sin^6>i3. 

Figure 3: Typical A4 Models. On the left (right), we plot sin^^ia as a function of 
sin^^i2 (sin^^23), following eq. (2.22). The dashed-black lines represent the 3a values for 
the mixing angles from the Fogli et ai fit [7]. Only the NH data sets is shown. The 
parameter ^ is taken equal to 0.076. The c^j^ parameters that multiply ^ are treated as 
random complex numbers with absolute values following a Gaussian distribution around 1 
with variance 0.5. 



The most general solution to this constraint can be parametrized in the following form: 

(X y — uj y + w \ 

y — w X + z + w y — z , (2.24) 
y + w y — z x + z — w J 

We see that w describes the deviation of m^^ from the form associated to the Tri-Bimaximal 
mixing, see eq. (2.12). The matrix rrij^ can be diagonalized in two steps. First we transform 
rUjy by a Tri-Bimaximal rotation: 

(X — y ^/Sw \ 

x + 2y . (2.25) 

V^w x-y + 2z J 

Second, we perform a unitary transformation in the (1,3) plane: 

fa r \ 

V = \ 10, |tt|2 + |^'|2 = 1^ (2.26) 

\ a* J 

V^m'^V = mf"^ . (2.27) 

The exact rotation is given by: 



2af uv*{u*-v) _ _2V3w , _ 2^3 



w 



with u = , and v = . (2.28) 



|q,|2 _ |^/|2 |y|2 _ |^|2 X — y ' X — y + 2z 
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The unitary matrix that diagonahzes rrij, is 



VV^a 1/73 VV^C 

UtbV= I -a/A/6 + C'7y2 1/73 -0^772-^776 
-a/76-C'772 1/73 +a772-CV76 



(2.29) 



Such a mixing pattern is very interesting because the observed 9is can be reproduced 
by choosing of order 0.1 and the predicted value of sin^ deviates from 1/3 only by 
terms of order To preserve these properties corrections from the charged lepton sector 
should be small compared to This can be realized by adopting two different expansion 
parameters <C for the charged lepton sector and for the neutrino sector. A model 
along these lines has been built in Ref. [61]. The setup is arranged in such a way that 
($e) and ($^) satisfy decoupled equations up to NLO so that it is possible to achieve 
($e) <C and to maintain the property of eq. (2.7) up to NLO. Moreover ($^) couples 
to charged leptons only at the NNLO so that the dominant source of corrections to the 
neutrino mixing pattern is 5m^J\ which in turns, being dominated by is invariant 

under S. 

In eq. (2.29) it is not restrictive to choose a real and positive and we have: 



Sep 
sin 9is 



arg 



^12 ^13 



(2.30) 
(2.31) 



sin^ e 



1 



12 



3-2|f|2 3 



l + lw\'-lneici, + ci,)i 



(2.32) 



sin^ 



1 + 



^/3a 



1 + 



^/3a 



23 



3a2 



1 



(2.33) 



where we have also included the effects coming from the diagonalization of the charged 
lepton sector as in eq. (2.19), to first order in ^. The second equality shows the result 
expanded in powers of |^'|, to the order In these models |^'| is chosen to be of order 
0.1, larger than ^ so that the contribution of eq. (2.19) can be neglected, and the lepton 
mixing matrix is very close to Utb V- It is interesting to note that if we neglect the 
corrections proportional to ^, we have exact correlations between the reactor angle and the 
other two angles^: 



sin^ e 



12 



3(1 - sin^ ^13) 



sin^ 



11-^ 

23 = - + Sm 013 COS dcp 



(2.34) 



^It has been shown in Ref. [69], from general group theoretical considerations, that these correlations 
are a general feature of flavour models when the symmetry group of the charged lepton (neutrino) mass 
matrix is Z3 {Z2). 
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The first expression shows that the unitary transformation V always increases the solar 
angle from the TB value, while the preferred la interval is below the TB prediction, see 
Fig. 1. This is a small effect, of second order in ^13, that can be compensated by the 
corrections proportional to ^. The second correlation involves the Dirac CP phase and is 
particularly interesting considering the recent hint of a CP phase close to tt for the NH 
case. In Fig. 4, we graphically compare this second expression with the present data for 
the NH case: when considering the la (2a) ranges for the mixing angles, one sees an 
indication that cos^^p lies in the interval [—1, —0.5], while no indication arises when the 
3a error band for sin^ 6^23 is taken into account. Although these results for the CP phase 
is modified by the inclusion of the subleading ^ contributions, these correlations will allow 
an interesting test for such models once Sep is measured and the precision on sin 623 is 
improved. 



Figure 4: Special A4 Models. Contour plot for 1, 2, 3a values of sin^ 623 in the parameter 
space sin^ Ois^cos 5cp according to the values for the NH in Table 1. The Yellow, Orange, 
Red regions refer to the data at la, 2a, 3a, respectively. 

In Fig. 5, we study the success rate to reproduce all the three mixing angles inside their 
corresponding 3a error ranges, as a function of |(^'|. The parameters have been chosen such 
that is a real number in [0.005, 0.06] and c^j are random complex numbers with absolute 
values following a Gaussian distribution around 1 with variance 0.5. The value of \^'\ that 
maximizes the success rate for both NH and IH is 0.184. The corresponding success rate 
is much larger in these models (~ 55%) than for the typical A4 models. 

We analyze quantitatively the deviations in eqs. (2.32) and (2.33) and their correlations 
in Fig. 6: in the plots on the left (right) column, we show the correlations in eqs. (2.32) 
and (2.33) between sin^ ^13 and sin^ 612 or sin^ ^23, respectively. The parameters have been 
chosen such that is a real number in [0.005, 0.06]; is a complex number with absolute 
values equal to 0.184; the parameters c^j are random complex numbers with absolute values 
following a Gaussian distribution around 1 with variance 0.5. In the plots we show only 
the NH case. The IH case is similar. For this choice of the parameters, the model can well 
describe all three angles inside the corresponding 3a interval, and its success rate is much 
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Figure 5: Special A4 Models. Success Rates as a function of the parameter \^'\. The 
parameters have been chosen such that is a real number in [0.005, 0.06] and c^j are 
random complex numbers with absolute values following a Gaussian distribution around 1 
with variance 0.5. In Cyan the NH and in Orange the IH. The value of |^'| that maximizes 
the success rate for both NH and IH is 0.184. 

larger than that of the typical TB models, as turns out by comparing Figs. 2 and 5. 



0.3 OJ 0.6 0.7 0^ 

(a) Correlation between sin^ O12 and sin^ 6^13. (b) Correlation between sin^ O23 and sin^ 6^13. 

Figure 6: Special A4 Models, sin^ 6^13 as a function of sin^ 6^12 (sin^ ^23) is plotted on 
the left (right), following eqs. (2.32) and (2.33). The dashed-black lines represent the 3a 
values for the mixing angles from the Fogli et ai fit [7]. Only the NH data sets is shown. 
The parameter is a real number in [0.005, 0.06]; is a complex number with absolute 
values equal to 0.184; the parameters c^j are random complex numbers with absolute values 
following a Gaussian distribution around 1 with variance 0.5. 




2.2 S4 Models 

In this section we refer to a SUSY model based on the fiavour symmetry Gf = S4 x 
Z4 X U{1) [52], but we keep the presentation slightly more general, to embrace a wider 
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class of possibilities [53-55]. The group S4 admits two generators S and T satisfying 



= {STf =T^ = 1. (2.35) 

Its irreducible representations are two singlets 1 and 1^, a doublet 2 and two triplets 3 and 
3^ In one of the two triplet representations of 6^4, S and T can be chosen as: 
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(2.36) 
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In the class of models considered here the electroweak SU(2) lepton doublets / transform as 
a triplet 3, the electroweak singlets e^, /x^ and transform as singlets and the electroweak 
Higgs doublets H^^ are invariant. If present, the right-handed neutrinos transform as a 
triplet 3. Similarly to the previous class of models, at the LO and in the SUSY limit the 
VEVs of $e and are determined by two decoupled equations and satisfy 

T'($e) = ($e) , S{^,) = [LO] . (2.37) 

The transformation = iT generates a Z4 subgroup of the flavour group Gf [52]. The 
charged lepton mass matrix rrie is given by 



nip = m 



+ (5m(^) + . . . (2.38) 



with the LO contribution mi^^ of the same type as the one considered before in eq. (2.9). 
In this case, r] represents a small parameter that breaks the flavour symmetry 6^4. At the 
LO the charged lepton mass matrix m^^\ only depending on (^e)? is diagonal and invariant 
under the action of the transformation T': 

r'^m(°)Wi°)r' = m(°)Wi°) . (2.39) 

We have = Z4, generated by T^ The hierarchical pattern <C 2/^ <C i/t is reproduced 
by operators of increasing dimensions contributing to y^-, and y^^ respectively. Also in 
this case we are including in mi^^ all terms arising from the LO (^e)? independently from 
the dimensionality of the operators that contribute to the charged lepton mass matrix. 

In the speciflc model of Ref. [52] a See-Saw mechanism produces a mass matrix for the 
light neutrinos m^^, given by: 

= m^^^ + Sm^J^ + . . . (2.40) 

where 




The parameters x, y and z are quadratic in the VEV of Hu and inversely proportional 
to the scale associated with the violation of the total lepton number. The LO term mi/^^ is 



14 



invariant under the x symmetry generated by S, eq. (2.37), and by the transformation 
^23 of eq. (2.13): 



m 



(0) 



(0) 



^23''^^ ^23 - • (2.41) 

Also in this case the Z2 symmetry represented by the matrix ^423 is an accidental symmetry. 
The matrix m>v^ of eq. (2.41) is the most general one invariant under both S and 7423. In 
the realization of Ref. [52] the parameters x, y and z are related by 



X 



r 

X 



(2.42) 



At the LO vfie^ is diagonal while vfiy^ is diagonalized by 



BM 



( 



\ 



\ 



I 



(2.43) 



This holds for any value of the parameters x, y and z since Ubm is the matrix that 
simultaneously diagonalizes S and A23. 

The departure from the bimaximal mixing depends on the subleading contributions 



bwbe , ^vfVu ' In the model of Ref. [52], such contributions are not generic. At the NLO 
the corrections to both the neutrino and the charged lepton sector are controlled by ($^) 
which preserves the LO alignment. The NLO correction to the charged lepton mass matrix 
is no longer invariant under T' and is of the following type: 



0(ye) 0{y:) 0{ye) 
0{yr) Oiyr) 
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where ^ is a small adimensional parameter given by the ratio between a flavon of the 
sector and A. The transformation needed to diagonalize mg is Vj'mJJe = mf*"^' where, to 
first order in ^ 
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-12 < 



^12 C ^13 ^ 
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(2.44) 



where c^j are complex number with absolute value of order one. In the neutrino sector 
after the inclusion of the NLO corrections the mass matrix has still the form of eq. (2.41) 
and is diagonalized by Ubm- The lepton mixing is U^Ubm and to first order in ^ we have 
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Figure 7: *S'4 Models. In (a), the success rate as a function of the parameter ^. The 
parameters 13 have been taken as random complex numbers with absolute value following 
a Gaussian distribution around 1 with variance 0.5. In Cyan the NH and in Orange the 
IH. The value of ^ that maximizes the success rate shown in the left plot is ^ = 0.172, for 
both the NH and IH. In (b), the success rate for the case c^g = as a function of ^ and 
cos Sep- The parameters c^2 been taken as a random complex number with absolute value 
following a Gaussian distribution around 1 with variance 0.5. The success rate reaches its 
maximum, ^ 20%, for cos^c'p = —1 and ^ = 0.18. Only the NH case is shown. The IH is 
similar. 



To properly correct the BM value of the solar angle to agree with the data, ^ is ex- 
pected to be 0{Xc)' Studying the success rate of having all three mixing angles inside 
the corresponding 3a ranges, we find that it is maximized for both the NH and IH when 
^ = 0.172, as shown in Fig. 7(a). In this case, the maximal success rate of about 2.6% 
is particularly small. The problem for this model is not to reproduce ^13 but rather to 
calibrate the correction to 6^12 for it to fall in its allowed window: we see from Fig. 7(a) 
that most of the scanning points spread out in a large interval of sin^ 612 between ^0.2 and 
^0.8. It may be interesting in this case to explore the possibility that one of the charged 
lepton mixing angles is dominant. For this to occur naturally an additional dynamical 
input would be needed. For the specific case c^3 = 0, we get a even more predictive cor- 
relation (still dependent, through 5c on the phase) among the solar and the reactor 
angle: 

sin^ 9i2 = ^ + sin 9is cos 5cp + 0(sin^ ^13) . (2.46) 

In this case, we study the dependence of the success rate from both ^ and the Dirac CP 
phase, in terms of co8 5cp' being dependent on two parameters, this success rate cannot 
be directly compared with the previous one, function of only ^; a trustworthy comparison 
requires to average among all the values of the success rate for a fixed ^. The corresponding 
plot for the NH case is shown in Fig. 7(b): the success rate reaches its maximum, ^ 20%, 
for ^ = 0.18 and co8 5cp = —1, for both the NH and IH cases. Only the NH case is shown 
in Fig. 7(b), while the IH is similar. The alternative possibility that = would lead to 
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similar results. 

It is interesting to investigate the origin of the preference for = —1 shown in 
Fig. 7(b). In Fig. 8, we graphically present the correlation in Eq. (2.46), neglecting the 
(^(sin^^ia) terms. The predicted value of sin^^i2 agrees with the experimental one only 
when C08 5c p is very close to —1. Although this correlation is modified once subleading 
contributions from both the neutrino and the charged lepton sectors are taken into account, 
still it will provide a strong test for such models. 

m \ \ ' " - 1 — I " T T T "I " ^ — M I 




OjOI 003 



Figure 8: S4 Models. Contour plot for 1, 2, 3a values of sin^ 612 in the parameter space 
sin^ ^13-cos 5c p according to the values for the NH in Table 1. The Orange and Red regions 
refers to the data at 2a and 3cr, respectively. No region corresponding to the data at la is 
present. 

We analyze quantitatively the expressions in eq. (2.45) and their correlation in Fig. 9(a), 
where 13 have been taken as random complex numbers with absolute value following a 
Gaussian distribution around 1 with variance 0.5, while = 0.172. In Fig. 9(b), we analyze 
the specific case c^g = 0, where ^ = 0.18 and cos^^p = —1. In Fig. 9, only the NH case is 
shown. The IH case is similar. Comparing Fig. 7(a) with Figs. 2 and 5, we can see that 
these 5^4 models are strongly disfavoured with respect to the A4 ones, and especially with 
respect to the special A4 models. 

In conclusion, from the point of view of reproducing the observed values of the mixing 
angles in a natural way, the A4 models of the Lin type provide a most efficient solution. In 
the next sections we will study the performance of the different models for LFV processes. 

3 Lepton Flavour Violation 

In this section we discuss the implications on LFV processes of the three classes of 
models described in the previous section. In particular, given the stringent upper bound 
on BR{ii 67), we focus on radiative lepton decays. By working in the super-CKM basis, 
where all kinetic terms are canonical and lepton mass matrices have been diagonalized 
through unitary transformations acting on the whole supermultiplets, the only sources of 
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(a) ci3 + (b) ci3 = 

Figure 9: S'4 Models, sin^ ^13 as a function of sin^ Q\2 is plotted, following eq. (2.45). The 
dashed-black lines represent the 3a values for the mixing angles from the Fogli et ai fit [7] . 
Only the NH data sets is shown. On the left, the parameters c^j are random complex 
numbers with absolute values following a Gaussian distribution around 1 with variance 
0.5, while ^ = 0.172. In the right, cl^ = 0, ^ = 0.18 and cos^cp = -1- 



LEV are the ofl[-diagonal terms of the slept on mass matrices. In our models these terms are 
much smaller than the corresponding diagonal entries and we can use the mass-insertion 
(MI) approximation [79^81] to illustrate the qualitative behavior of the model predictions. 
The quantitative results shown in our plots have been obtained by using complete one-loop 
results, which can be found for example in Refs. [82-88]. The normalized branching ratios 
Rij for the LFV transitions k Ij^ 

BR{k ^ IjViVj) ^ ' 

can be written as 

R^^=T^^m'^\m. (3.2) 

^ Fill SUSY 

At the LO in the MI approximation, the amplitudes A^l and are given by: 

J^' (3.3) 

— ^RRy^ij RR + ^LR [^ij LR 

rrii 

where rrii are the charged fermion masses and ace {C^ C = L, R) are dimensionless func- 
tions of the SUSY parameters msusY-> ^1,2, tan /5, renormalized at the electroweak scale. 
Here a common value in the diagonal entries of both LL and RR blocks of the slepton mass 
matrices has been assumed at the electroweak scale and rrisusY denotes the average mass. 
Such assumption, often made at the cut-ofl[ scale in constrained versions of the MSSM, is 
spoiled by running effects and does not hold any more at the electroweak scale. We will 
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discuss this effect in the next subsection. The simphfied framework considered here is suf- 
ficient to correctly describe the relation between Rij and the expansion parameters ^ and 
In our conventions the explicit expression of ace is given in appendix A. Their typical 
size is one tenth of 5'^/(167r^), g being the SU(2)i, gauge coupling constant. To appreciate 
the relative weights of the contributions in eq. (3.3), we list in Tab. 3 the expressions and 
the numerical values of the functions acc^ in the limit /x = Mi^2 = "^susy- As one can see, 
in this limit the dominant coeflScient is a^^, which is larger than a^L = cllr ^ factor 
7 -i- 54, and larger than qrr by a factor 4 -i- 14, depending on tan /3 = 2 ^ 15. This range 
of tan /3 is taken as representative for the models under consideration. More precisely, the 
parameter tan /3 is related to the expansion parameter ry, the mass of the r lepton and the 
T Yukawa coupling by: 

tan/5^^E^, (3.4) 

where Vew ^ 246 GeV is the EW Higgs VEV. For typical A4 models the parameter r] ^ ^ 
is of order 0.1. For special AA model, rj ^ ^ is smaller than and we will use the range 
0.007 ^ ^ ^ 0.05. For the S4 models, we have 77 0.08 to correctly fit the charged lepton 
masses. Requiring to be of order one, 1/3 ^ ^ 3, we have the following allowed 
ranges for tan /3: 

3 < tan /5 < 30 Typical 

2 < tan;5 < 15 Special A^ (3.5) 

3 < tan/3 < 24 ^4 



ClLL 


^ ^,[1 3(l + 4c)tan^6'vF + 4f4 + 5tan^6'vF') tan/?] 
240 IGtt^ L V / J 


+(2.0^ 16.3) 




1 , 2/3 

r tan vw 

12 W 


0.30 


dRR 


1 ^2 

60 167r2 ^ ^""""^^ 


-(0.5^ 1.2) 



Table 3: Coefficients ace characterizing the transition amplitudes for /x ^ 67, r ^ 67 
and T fij^ in the MI approximation and by taking the limit the /i = Mi^2 = '^susy- 
Numerical values are given in units of 5'^/(1927r^) and using sin^^^^ = 0.23, c = 1, and 
tan/3 = 2 ^ 15. The parameter c is a model dependent quantity of order one. 



Finally, {5ij)cc' parametrize the Mis and are defined as: 

[Oij)cc' - —2 , l^-oj 

where the different blocks of the slepton mass matrices mlcc' evaluated in the super- 
CKM basis, denoted by the hat. There are two main types of contributions to the Mis 
{5ij)cc''' the first one comes from local operators, bilinear in the slepton fields, with in- 
sertions of the fiavons $e,z/- They are invariant under the fiavour symmetry and represent 
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the counterpart, in the sparticle sector, of the operators that generate lepton masses in the 
theory. After the breaking of the flavour symmetry, we get slepton mass matrices expanded 
up to a certain order in ($e^^)/A, depending on the highest dimensionahty of the operators 
included: 

= mg + Srhl + Srhl + ... (3.7) 

with 5ml order (($e^^)/A)^. 

The second contribution to the Mis in eq. (3.6) comes from the renormahzation group 
evolution (RGE). When neutrino masses are generated by the See-Saw mechanism, the 
evolution spans two main regions. The flrst one goes from the cut-ofl[ scale A where slepton 
masses are generated'', down to the lightest right-handed neutrino mass Mi. In our models 
the mass scale associated to right-handed neutrinos coincides with the flavour symmetry 
breaking scale ($^^) and, in a leading logarithmic approximation, we get a contribution to 
the Mis proportional to y'^/ (167r^) x log(($^^)/A), y representing a typical neutrino Yukawa 
coupling. Below Mi the right-handed neutrinos decouple and the running is only aflFected 
by the lightest degrees of freedom, those of the MSSM. This part of the RGE occurs 
whether or not neutrino masses originate from the See-Saw mechanism and gives rise to 
contributions to the Mis proportional to y^/(167r^) x log{{^jy) /msusY)^ Vt denoting the r 
Yukawa coupling. 

We now discuss in turn the different contributions to the Mis for the three classes of 
models. 

3.1 Mass Insertions from Local Operators 

Local operators giving rise to slepton masses can be constructed with standard tech- 
niques, see Refs. [89-100]. The inclusion of non-renormalizable operators with insertions 
of the flavon flelds $e^^ also affects kinetic terms both in the fermion and sfermion sector, 
providing an additional source of flavour violation. Here we list the slepton mass matrices 
in the super-CKM basis where kinetic terms have been set in the canonical form by means 
of appropriate transformations and where fermion mass matrices have been made diago- 
nal through unitary transformations acting on the whole supermultiplet. These slepton 
masses refer to speciflc models, taken as representative of the classes analyzed above. For 
typical A4 models we will refer to the construction in Ref. [71], for special A4 models we 
will consider the model of Ref. [61] and flnally S4 models are exemplifled by the model 
in Ref. [52]. The results given below have been obtained under the assumption that the 
underlying parameters are real. 

^Depending on the specific type of SUSY breaking mechanism slepton mass generation can occur at a 
scale smaller than A. Here we assume that slepton masses are produced at the highest scale, like in gravity 
mediated SUSY breaking scenarios. 
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3.1.1 Typical A4 Models 

For typical A4 models we get 
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where dots stand for negligible SUSY contributions. In m^^ we have neglected a con- 
tribution which arises if the F component of the flavon supermultiplets acquires a VEV. 
Such a contribution depends on the SUSY breaking mechanism and vanishes under mild 
assumptions [96, 101, 102]. A similar contribution will be neglected also in the special A4 
and in the S4 models discussed below. We have 
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and 



A^i = a^LOie) + an^Oie) , 

Al = anR'^Oie + anJ^Oie 

rrij rrij 



(3.12) 



From 0.1 and the numerical values of ace we see that the amplitude A^l is the 
dominant one. We approximately have 



J-Hj — 



^F^SUSY 



(3.13) 



and we expect the branching ratios of the three transitions to be of the same order of 
magnitude: 



at variance with the predictions of most of the other models, where, for instance, i?^e/i?^^ 
can be much smaller than one [103-106]. 
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3.1.2 Special A4 Models 

For the special A4 model in Ref. [61] we get 
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Neglecting the subdominant contribution from the yl^ amphtude, we have 
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and we expect the branching ratios of the three transitions to be of the 
magnitude: 



3.1.3 ^4 Models 

For 5^4 models we get: 
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We have 
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Neglecting the subdominant contribution from the ^4^ amphtude, we have 
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with a suppression of the rate t ^ ji^ relative to /i ^ 67 and r ^ 67 by a factor of 

Rre . (3.28) 



Before a more quantitative illustration of these results, we discuss the effects induced by 
the RGE from the scale of flavour symmetry breaking down to the electroweak scale. 

3.2 Mass Insertions from low-energy RGE 

This set of corrections is common to all models irrespective of the assumed existence 
of RH neutrinos, i.e. both with or without See-Saw. The running of the slepton mass 
parameters from the scale ($e,i/) down to the scale msusY can be estimated in a leading 
logarithmic approximation. The largest effect is a correction to the matrices rh\^ and 
coming from electroweak gauge interactions and is proportional to the identity matrix in 
flavour space. One flnds that the diagonal elements increase in the running from the cutoff 
scale down to the electroweak scale. This effect is taken into account in our numerical 
study. As for the off diagonal entries the largest corrections are proportional to the square 
of the T Yukawa coupling. Since ^^/(IGtt^) 3 x (10"^ ^ 10"^) for tan^S = (2 ^ 15), even 
in the presence of the large factor \og{{^ ^.u) I'^susy) ^ 30, these corrections are negligibly 
small compared to the contribution from the local operators discussed above. We can 
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conclude that the corrections to the off-diagonal entries of the soft mass matrices induced 
by the RGE from ($e,z/) down to rrisusY are either negligible or could be absorbed in the 
parametrization given in the previous section. 

The present bound on the branching ratio of /x ^ 67 from the MEG collaboration [64] , 
BR[ii ej) < 2.4 X 10~^^, leads to strong constraints on the parameter space of the 
models that we are considering in this work. The supersymmetric parameters that are not 
constrained by the flavour symmetry, such as the soft SUSY mass scales and the gaugino 
and Higgs(ino) sectors, are flxed by our choice of a SUGRA framework: mo and M1/2 are 
the common masses of scalar particles and gauginos at the GUT scale. Thus, at the scale 
A = 2 X 10^^ GeV, 

Mi(A) = M2(A) = Mi/2, (3.29) 

where Mi are the SU{2) x U{1) gaugino masses. The effects of the RG running lead at low 
energies to the following masses for the gauginos 

M^{m^) - ^^M,(A) M,(m^) ^ ^^M,(A) , (3.30) 

where ai = gf/^i^ {i = 1,2) and ai{A) = a2{A) 1/25. We have seen that, among 
the RG running effects on the soft mass terms, only those from the electroweak gauge 
interactions are relevant. At the cut-off scale the LL and RR blocks of the slepton mass 
matrices are given by the eqs. (3.8)-(3.10), (3.15)-(3.17) and (3.22)-(3.24) with tususy 
identifled with mo. The RGE due to electroweak gauge interactions leaves the off-diagonal 
entries essentially unaffected, while the diagonal elements at the weak scale, denoted by 
Tn\j^{mw)^ are given by: 

ml{mw) ^m^ + 0.54Mi%, 
Notice that this effect modifles the previous estimates of the LL{RR) Mis by a factor 

2 

1/2 



ml/ {ml + 0.54(0.15)M^/2)5 thus providing an additional suppression when M1/2 is larger 
than mo- 

Furthermore, the parameter /i is flxed through the requirement of correct electroweak 
symmetry breaking^ 



^ _mi l + 0.5tan^/j , 0.5 + 3.5 tan^ /j 
1^1 - - ^'^o tan2/3-l tan2/5-l ' 



SO that fi is determined by mo, M1/2 and tan^ up to its sign. We recall that in our model 
the low energy parameter tan /5 is not a free parameter, as shown in Eq. (3.4). We compare 
the models at tan/3 = 2 and tan/3 = 15. In our numerical analysis, we have assumed 



^The general definition of the parameter /i is 



that reduces to the expression in Eq. (3.33) once considering that in the SUGRA framework the soft Higgs 
mass parameters are also given by mg at the high energy scale, m|^^(A) = m'jj^{A) = m? 



0- 
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that the parameters on the diagonal of the slepton mass matrices (m^^ {rnlpijijK 
are positive in order to get positive definite square-masses and to avoid electric-charge 
breaking minima and further sources of electroweak symmetry breaking. The absolute 
value of the 0{1) parameters is varied between 1/2 and 2. Furthermore, we have imposed 
the conditions that the lightest chargino has a mass larger than 100 GeV and that the 
lightest neutralino is the lightest supersymmetric a particle (LSP). 

The results for the typical models are illustrated in Fig. 10, where the parameter 
^ is taken equal to 0.076, consistently with the analysis in Sec. 2.1.1. For tan/3 = 2 and 
mo = 200 GeV, Fig. 10(a), almost aU the points are excluded for M1/2 < 400 GeV. The 
corresponding supersymmetric spectrum is rather light: for M1/2 = 400 GeV, the lightest 
neutralino has a mass of ^ 156 GeV, the lightest chargino of ^ 306 GeV, the lightest 
LH charged slepton is in the range [230, 500] GeV and the lightest RH charged slepton 
in the range [160, 350] GeV. These mass values are not yet excluded by the LHC (the 
limits in the EW sector are not very strong) and they only mark the lower edge of the 
region allowed by the fi ^ ej bounds. Increasing mo up to 5000 GeV, Fig. 10(b), we 
see that the MEG bound is well satisfied in this model for the whole plotted range of 
Ml/2. The corresponding supersymmetric spectrum is heavier: again for M1/2 = 400 GeV, 
while the lightest neutralino and chargino have masses very similar to the previous case, 
^ 158 GeV and ^ 315 GeV, respectively, the lightest LH and RH charged slepton masses 
are much higher, being in the range [3850, 7070] GeV and [1800, 6260] GeV, respectively. 
Increasing mo from 200 GeV up to 5000 GeV does not correspond to a uniform decrease 
of the branching ratio, as a function of M1/2. The plot with mo = 5000 GeV is much 
fiatter compared to the plot where mo = 200 GeV. This is due to the approximate factor 
ml/ {ml + 0.54M^/2) entering the dominant LL mass insertion, which, for the values of 
Ml/2 used in our plots, is sharply decreasing for mo = 200 GeV while it is slowly varying 
and close to one for mo = 5000 GeV. This is a general feature reproduced by all the models 
considered here. 

Increasing the value of tan^S, Figs. 10(c)-(d), the number of points below the MEG 
bound decreases, but the previous considerations are approximatively still valid: in partic- 
ular notice that for mo = 5000 GeV, there are always points satisfying the MEG bound, 
even if the largest number of them falls in the excluded region, especially for smaller M1/2. 
Since the dominant contribution to the branching ratio comes from the LL mass insertion 
which is proportional to tan the branching ratio is to a good approximation proportional 
to tan^ /3, as we can see from the plots. Also this feature is common to all models. 

For special A4 models the parameter is taken equal to 0.184, as analyzed in Sec. 2.1.2. 
The results of our numerical study considering the constraint from BR{ii 67) for the 
mo — Ml/2 parameter space are shown in Fig. 11. The plots in Fig. 11 are very similar to 
those in Fig. 10. Indeed also for this case, for tan/3 = 2 and mo = 200 GeV, Fig. 11(a), 
almost all the points are excluded for M1/2 ^ 400 GeV. The supersymmetric spectrum 
is also quite similar: the only differences are in the range of masses that the lightest LH 
and RH charged sleptons can span, [260, 500] GeV and [180, 340] GeV, respectively. By 
increasing mo up to 5000 GeV, Fig. 10(b), we can see that for the whole range of M1/2, 
the MEG bound is well satisfied in this model, while the corresponding supersymmetric 
spectrum is heavier. Special and typical A4 models have a dominant LL mass insertion 
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A/,y2 (GeV) Mm2 (GeV) 



(a) tan/? = 2 and mo = 200 GeV (b) tan 13 = 2 and mo = 5000 GeV 




A/,/2 (GeV) Ml/2 (GeV) 

(c) tan/3 = 15 and mo = 200 GeV (d) tan/3 15 and mo 5000 GeV 

Figure 10: Typical A4 Models. Scatter plots of BR{fi 67) as a function of M1/2, for 
different values of tan /3, and rriQ. The parameter is chosen as 0.076 in order to maximize 
the success rate of this model. The horizontal line shows the current MEG bound. For 
Blue (Red) points the LSP is the lightest neutralino (stau). The percentage in each plot 
refers to the number of Blue points that satisfy the MEG bound over the total number of 
points. 

proportional to (^'^ and (^^, respectively which represents the main difference between the 
two models, as far as radiative charged lepton decays are considered. The optimal values 
of ^ and l^'l differ by a factor of about two and we expect that the branching ratios of the 
two models should differ by about one order of magnitude, for fixed values of the other 
parameters. This effect is barely visible in our plots, due to the spread of the predictions 
caused by the unknown order-one coefficients. 

Increasing the value of tan /5, the number of points below the MEG bound decreases, but 
the previous consideration are approximatively still valid: the most interesting difference 
is in Fig. 11(c), for mo = 200 GeV, where almost all the points with M1/2 < 1000 GeV are 
excluded; furthermore, in Fig. 11(d), for mo = 5000 GeV, there are always points satisfying 
the MEG bound, even if the largest part are excluded, especially for smaller M1/2. 

Finally, in the S4 model the parameter ^ is taken equal to 0.172, as required to maximize 
the success rate of these models to arrange the three mixing angles in the corresponding 
3a ranges. The results are displayed in Fig. 12. The plots in Fig. 12 are very similar to 
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A/,y2 (GeV) Af 1/2 (GeV) 



(a) tan/? = 2 and mo = 200 GeV (b) tan^S = 2 and mo = 5000 GeV 




A/,y2 (GeV) Ml/2 (GeV) 

(c) tan/3 = 15 and mo = 200 GeV (d) tan/3 15 and mo 5000 GeV 

Figure 11: Special A4 Models. Scatter plots of BR{fi 67) as a function of M1/2, for 
different values of tan /3, and rriQ. The parameter |^'| is chosen as 0.184 in order to maximize 
the success rate of this model. The horizontal line shows the current MEG bound. For 
Blue (Red) points the LSP is the lightest neutralino (stau). The percentage in each plot 
refers to the number of Blue points that satisfy the MEG bound over the total number of 
points. 

those in Fig. 11 and the same comments also apply here. More interestingly, in all the 
plots the number of points satisfying the MEG bound is much smaller, especially for large 
mo. In particular for tan /5 = 15 and mo = 5000 GeV, only the ^ 10% of points correspond 
to a BR smaller than the MEG bound. The relatively larger branching ratio for /x ^ 67 
predicted by the S4 model is also a consequence of the scaling of {5i^e)LL with respect to ^\ 
such a scaling is linear in 5^4, while it is quadratic in the typical models. Moreover the 
optimal value of ^ in is larger than in typical models. This explains the enhancement 
by two order of magnitude of the 6^4 prediction compared to the typical one. 

3.3 Mass Insertion from high-energy RGE 

When neutrino masses originate from a Type I See-Saw mechanism, there is an extra 
contribution to the running of the slepton mass matrices, originating from loop diagrams 
with the exchange of right-handed neutrinos. Such a contribution is only relevant in the 
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A/,y2 (GeV) Mm2 (GeV) 

(a) tan/? = 2 and mo = 200 GeV (b) tan^S = 2 and mo = 5000 GeV 




A/,y2 (GeV) Ml/2 (GeV) 

(c) tan/3 = 15 and mo = 200 GeV (d) tan/3 15 and mo 5000 GeV 

Figure 12: 84^ Models. Scatter plots of BR{ii ej) as a function of M1/2, for different 
values of tan/3, and mg. The parameter ^ is chosen as 0.172 in order to maximize the 
success rate of this model. The horizontal line shows the current MEG bound. For Blue 
(Red) points the LSP is the lightest neutralino (stau). The percentage in each plot refers 
to the number of Blue points that satisfy the MEG bound over the total number of points. 

energy range from the cut-off A of the theory down to the right-handed neutrino masses 
Mk (k = 1,2,3). By focussing on the LL block m^^^ of the slepton mass matrix, whose 
non-diagonal entries dominate the rates of the processes under consideration, in the leading 
log approximation we have [107-109] 

(m,2^^)^^.^-^(3m^ + A^)5](y;%log(A) (i;),^., (3.34) 

where is the SUSY breaking parameter characterizing the size of the trilinear scalar 
mass term for sleptons and the matrix denotes the neutrino Yukawa couplings in the 
basis where the mass matrices for charged leptons and right-handed neutrinos have been 
diagonalized. The above expression for ml^^ holds at the scale equal to the lightest right- 
handed neutrino mass. Below that scale, right-handed neutrinos do not affect the running 
any more. In the models we are considering neutrino Yukawa couplings are of order one 
and even for relatively small ratios A/Mk ^ 100, we may easily get contributions to the 
LL mass insertions of order 0.1. 
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It is interesting to note that if neutrinos transform as an irreducible triplet of the flavour 
symmetry then, at the LO in the flavour symmetry breaking parameters, the right-hand 
side of eq. (3.34) can be expressed in term of the light neutrino masses and the matrix 
elements of the lepton mixing matrix U. We have 



% = kU^ + ... (3.35) 



and 



2 



Ui2 log — + Uis log — [/-g 
mi mi 



+ ..., (3.36) 



where is a constant of order one and dots stand for non-leading contributions, suppressed 
by powers of ($)/A. This result is completely general. It applies to any model with a flavour 
symmetry, independently of the speciflc flavour group G/, provided neutrinos are assigned 
to an irreducible triplet of the group, like in the models under consideration. First note 
that the invariance of the theory under transformations of implies 

p{g) YX P{9)^ = yX , (3-37) 

p denoting the irreducible triplet representation under which the left-handed leptons trans- 
form. The combination Y^Y^, commutes with each group element p{g) and, by the Shur's 
First Lemma, Y^Y^ either vanishes, a case that we exclude, or is proportional to the unit 
matrix. We conclude that Y^ is proportional to a unitary matrix. This result holds in 
any basis, since the invariance of the theory is a basis-independent property. If now we go 
to the basis where charged leptons and right-handed neutrinos are mass eigenstates, the 
See-Saw relation reads: ^ 

m, = ^yjM-^i; (3.38) 

and we recognize that in this basis the unitary matrix to which Yy is proportional should 
coincide with [/^ We obtain eq. (3.35) and 

= j^,{m^Ua, (3.39) 

By making use of eqs. (3.35), (3.39) and (3.34) we immediately get the result in eq. (3.36). 
Thus, up to sub-leading corrections and up to the unknown order-one parameter A:, the LL 
mass insertions are completely determined by neutrino masses and mixing parameters. At 
the LO the Mis do not depend on the cut-ofl[ scale A, but only on the ratios between light 
neutrino masses. For a degenerate neutrino spectrum the LO Mis vanish. On the contrary 
the largest Mis are obtained when the spectrum is hierarchical. 

We have considered the type I See-Saw version of the previous models and we have 
evaluated the slepton mass matrices by including the eflFects on the running due to right- 
handed neutrinos. Our plots and our numerical results have been worked out at the NLO 
in the symmetry breaking parameters. RGE equations are solved numerically, making use 
of full one-loop beta functions. The evolution starts at A = 2 x 10^^ GeV, where we assume 
the pattern dictated by the eqs. (3.8)-(3.10), (3.15)-(3.17) and (3.22)-(3.24) with tususy 
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M,j2 (GeV) Mm2 (GeV) 



(a) tan/? = 2 and mo = 200 GeV (b) tan 13 = 2 and mo = 5000 GeV 




A/,/2 (GeV) Ml/2 (GeV) 

(c) tan/3 = 15 and mo = 200 GeV (d) tan/3 15 and mo 5000 GeV 

Figure 13: Typical A4 Models. Scatter plots of BR^ji ej) with running effects due to 
right-handed neutrinos, as a function of M1/2, for different values of tan/3, and mo. Light 
neutrinos have normal mass ordering witn mi = 4.4 meV. The parameter ^ is chosen as 
0.076 in order to maximize the success rate of this model. The horizontal line shows the 
current MEG bound. For Blue (Red) points the LSP is the lightest neutralino (stau). The 
percentage in each plot refers to the number of Blue points that satisfy the MEG bound 
over the total number of points. 

identified with ttiq. Going down to the electroweak scale the off-diagonal entries of the 
slepton mass matrices are modified by the running due to both right-handed neutrinos and 
by the MSSM degrees of freedom. We display our results for the case of normal ordering. 
Similar considerations hold when the neutrino mass ordering is inverted. For the lightest 
neutrino mass we chose the smallest value allowed by the models under consideration: 4.4 
meV for the typical A4 models, 0.4 meV for the special A4 models and 1 meV for S4. These 
values can be estimated by analyzing the neutrino masses in the See-Saw version of the 
models in Refs. [71], [61] and [52], considered in the previous subsections. Choosing the 
smallest value of mi enhances the ratios 1712/1711 and ms/rrii in eq. (3.36) and maximizes 
the effect of the running due to right-handed neutrinos. 

In Figs. 13, 14 and 15 we plot the branching ratio for // ^ 67 as a function of M1/2 
for the same values of mg and tan^S shown in Figs. 10, 11 and 12 to allow for a direct 
comparison between the two cases, with and without right-handed neutrinos. As a general 
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M,i2 (GeV) Mm2 (GeV) 



(a) tan/? = 2 and mo = 200 GeV (b) tan 13 = 2 and mo = 5000 GeV 




A/,y2 (GeV) Ml/2 (GeV) 

(c) tan/3 = 15 and mo = 200 GeV (d) tan/3 15 and mo 5000 GeV 

Figure 14: Special A4 Models. Scatter plots of BR{ii ej) with running effects due to 
right-handed neutrinos, as a function of M1/2, for different values of tan/3, and mg. Light 
neutrinos have normal mass ordering witn mi = 0.4 meV. The parameter \^^\ is chosen as 
0.184 in order to maximize the success rate of this model. The horizontal line shows the 
current MEG bound. For Blue (Red) points the LSP is the lightest neutralino (stau). The 
percentage in each plot refers to the number of Blue points that satisfy the MEG bound 
over the total number of points. 

trend, the inclusion of the running due to right-handed neutrinos enhances the branching 
ratio by a factor between about one and two orders of magnitude. The allowed region 
in the parameter space of the models shrinks, as indicated by the percentage of points 
that satisfy the current MEG bound. The largest effect occurs for the special A4 model, 
also due to the small value mi = 0.4 meV, which, for normal hierarchy, enhances the 
contribution in eq. (3.36). In the cases of typical A4 models and S4 models, the effect of 
right-handed neutrinos is similar. Notice that the slope of the plotted regions for mo = 200 
GeV is much steeper than for mo = 5000 GeV which increases the impact of right-handed 
neutrinos in the latter case. This is particularly visible in the case of special A4 models. If 
we compare the panel (d) of Figs. 11 and 14, where mo = 5000 GeV and tan^S = 15, we 
see that the points are lying on an almost horizontal strip close to the MEG bound and 
the enhancement due to right-handed neutrinos is sufficient to displace almost all the band 
above the bound, thus excluding most of the points. In conclusion, running effects due 
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M,i2 (GeV) 

(a) tan/? = 2 and mo = 200 GeV 




A/,y2 (GeV) 

(c) tan/? = 15 and mo = 200 GeV 




Mm2 (GeV) 

(b) tan ;9 = 2 and mo = 5000 GeV 




Ml/2 (GeV) 



(d) tan/3 = 15 and mo = 5000 GeV 



Figure 15: *S'4 Models. Scatter plots of BR^ji ej) with running effects due to right- 
handed neutrinos, as a function of M1/2, for different values of tan/3, and mo- Light 
neutrinos have normal mass ordering witn mi = 1 meV. The parameter ^ is chosen as 
0.172 in order to maximize the success rate of this model. The horizontal line shows the 
current MEG bound. For Blue (Red) points the LSP is the lightest neutralino (stau). The 
percentage in each plot refers to the number of Blue points that satisfy the MEG bound 
over the total number of points. 

to right-handed neutrinos can significantly contribute to the off-diagonal terms of slepton 
mass matrices and can even dominate Mis, especially for a pronounced hierarchy in the 
light neutrino mass spectrum. There is a general reduction in the parameter space of the 
model and in particular for the special A4 model with the most hierarchical spectrum. A 
milder impact is expected for a neutrino mass spectrum close to the degenerate case. 

3.4 Correlation with the Muon g — 2 

The value found for the anomalous magnetic moment of the muon [110] 

a™ = 116592080(63) x 10"^^ (3.40) 

shows a 3.4 a deviation 

Sa^ = af^^ - a^^ = +302(88) x 10"^^ (3.41) 
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from the value expected in the SM [111-113] 

= 116591778(61) x 10"^^ . (3.42) 

It is an interesting question whether the presence of SUSY particles can account for this 
deviation, once the constrains from the branching ratio of the /x ^ 67 decay are taken into 
consideration. 

Following Refs. [114-117], we study the correlation between and BR^ji 67), for 
tan/5 G [2, 15] and mo, M1/2 G [200, 5000] GeV, while all the other parameters are treated 
according to the previous section. Only points corresponding to scenarios with the lightest 
neutralino being the LSP are shown. 




Figure 16: Correlation plots between 5a^ and BR{iJi 67). The value of tan/5 is taken in 
the range [2, 15], while mo, M1/2 are chosen between 200 and 5000 GeV. The values of ^ 
or \^'\ are those that maximize the success rates for the three models for the NH case. The 
colour of the points refer to the values of tan^S: 2 < tan/5 ^ 7 in blue, 7 < tan^S < 11 in 
red, 11 < tan/5 < 15 in green. The horizontal line corresponds to the MEG bound, while 
the vertical lines correspond to the measurements on 5a ^ at 3(7. 

As we can see from Fig. 16, in the whole parameter space of the three models it is not 
natural to reproduce the observed deviation of the muon anomalous magnetic moment, 
once we consider the present MEG bound on BR{fi 67). This is not a surprise, because 
the explanation of the 3.4a discrepancy needs small values of mo and M1/2 and larger 
values of tan/5, which, however, enhance the branching ratio of the radiative LFV decays. 



4 Conclusion 

The recent rather precise measurements of ^13 make our present knowledge of the neu- 
trino mixing matrix, except for the CP violating phases, sufficiently complete to consid- 
erably restrict the class of models that can reproduce the data. In spite of this progress, 
the range of possibilities for flavour models remains unfortunately quite wide. On the one 
extreme, the rather large value measured for ^13, close to the old CHOOZ bound, has 
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validated the prediction of models based on anarchy [118, 119], i.e. no symmetry in the 
leptonic sector, only chance, so that this possibility remains valid, as discussed, for exam- 
ple, in Ref. [120]. Anarchy can be formulated in a SU (5) ® U (1) context by taking different 
Froggatt-Nielsen [121] charges only for the SU (5) tenplets (for example 10 ^ (3, 2, 0), where 
3 is the charge of the first generation, 2 of the second, zero of the third) while no charge 
differences appear in the 5 (5 ^ (0,0,0)). Anarchy can be mitigated by assuming that it 
only holds in the 2-3 sector with the advantage that the first generation masses and the 
angle ^13 are naturally small (see also the recent revisiting in Ref. [122]). In models with 
See-Saw, one can also play with the charges for the right-handed SU(5) singlet neutrinos. 
If, for example, one takes 1 ^ (1, —1, 0), together with 5 ^ (2, 0, 0), it is possible to get a 
normal hierarchy model with ^13 small and also with r = Am^^^^^/Am^^^ naturally small 
(see, for example, Ref. [123]). In summary, anarchy and its variants, all based on chance, 
offer a rather economical class of models that are among those encouraged by the new ^13 
result. On the other extreme, stimulated by the fact that the data suggest some special 
mixing patterns as good first approximations (TB or BM, for example), models based on 
discrete fiavour symmetries, like or 6^4, have been proposed and widely studied. In 
these models the starting LO approximation is completely fixed (no chance), but the NLO 
corrections introduce a number of undetermined parameters. The recent data on ^13 and 
the MEG new upper bound on the LFV process fi ^ ej impose a reappraisal of these 
models, which we have attempted in this paper. In particular, the relatively large value of 
9is introduces a marked departure from the TB limit, while the values of 612 and 623 are 
very close to it. The challenge is to produce in a natural way a relatively large correction 
to 9is without affecting too much the other mixing angles. But one must pay attention 
that these larger corrective terms introduced to shift ^13 from the TB value could appear 
in the non-diagonal elements of the charged lepton (and s-lepton) mass matrix and could 
induce a too large fi ^ ej branching ratio. 

As a result of our analysis we find that, for reproducing the mixing angles, the Lin type 
A4 models have the best performance, as expected, followed by the typical A4 models, 
while the BM mixing models lead to an inferior score, as they most often fail to reproduce 
9i2. In the latter case, if only one complex parameter perturbs the BM pattern, there is 
a strict correlation among the solar and the reactor angle and the success rate increases 
considerably by selecting a CP violating Dirac phase close to tt. 

As for LFV processes we have addressed the problem by adopting the simple CMSSM 
framework. While this overconstrained version of supersymmetry is rather marginal after 
the results of the LHC searches, more so if the Higgs mass really is around mn = 125 
GeV, we still believe it can be used here for our indicative purposes. We find that the most 
constrained versions are the models that start with BM mixing at the LO because, in this 
case, relatively large corrections directly appear in the off-diagonal terms of the charged 
lepton mass matrix. The typical A4 models turn out to be the best suited to satisfy the 
MEG experimental bound, as the non-diagonal charged lepton matrix elements needed to 
reproduce the mixing angles are quite smaller. An intermediate, still rather good, score is 
achieved by the models of the Lin type, where the main corrections to the mixing angles 
arise from the neutrino sector. When the fit to the mixing angles and the bounds on LFV 
processes are combined, the A4 models emerge well from our analysis and in particular 
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those of the Lin type perhaps appear as the most reahstic approach to the data among 
the models based on discrete flavour groups that we have studied. As for the regions of 
the CMSSM parameter space that are indicated by our analysis the preference is for small 
tan /3 and large SUSY masses (at least one out of mo and mi/2 must be above 1 TeV). As a 
consequence it appears impossible, at least within the CMSSM model, to satisfy the MEG 
bound and, at the same time, to reproduce the muon g — 2 discrepancy. 
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A Expression of acc^ 

In this section we list the explicit expression of ace in our conventions: 
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(A.l) 

where ai,2 = Ml^/ml^sYi ^ = n'^/ml^sY and fi(c,n){x,y) = fi{c,n){x) - fi{c,n){y)- The 
functions fin{x) and fic{x), shghtly different from those in Ref. [124], are given by: 



finix) = 

f2nix) = 
hnix) = 

hc{x) = 

f2c{x) = 



-17x^ + 9x'^ + 9x-l + Qx^{x + 3) loga;)/(24(l - xf) 
-5x^ + 4x + 1 + 2x{x + 2) log x) /(4(1 - x)'^) 
l + 9x-9x^ -x^ + 6x{x + 1) loga;)/(2(l - xf) 
-x^ -9x'^ + 9X + 1 + 6x{x + 1) loga;)/(6(l - xf) 
-x^-4:X + 5 + 2{2x + 1) log x) /(2(1 - x^) . 



(A.2) 
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